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^ '. Abstract 

\ Based on the contractor renormalization group (CORE) method and the density matrix renor- 

maUzation group (DMRG) method, a new computational scheme, which is called the block density 

^ ■ matrix renormalization group with effective interactions (BDMRG-EI), is proposed to deal with 

c/3 , the numerical computation of quantum correlated systems. Different from the convential CORE 

^ • 

^ ' method in the ways of calculating the blocks and the fragments, where the DMRG method instead 

S ' 

I ' of the exact diagonalization is employed in BDMRG-EI, DMRG-EI makes the calculations of larger 

^ , blocks and fragments applicable. Integrating DMRG's advantage of high accuracy and CORE's 

^ ' advantage of low computational costs, BDMRG-EI can be widely used for the theoretical calcula- 

tions of the ground state and low-lying excited states of large systems with simple or complicated 
■ connectivity. Test calculations on a 240 site one-dimensional chain and a double-layer polyacene 

\ oligomer containing 48 hexagons demonstrate the efficiency and potentiality of the method. 
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I. INTRODUCTION 

The real space renormalization group (RSRG) method firstly proposed by Wilson in 
1975, is a variational scheme that can solve large systems through truncating the Hilbert 
space successively. Despite RSRG's great success in weakly correlated systems, lots of the- 
oretical calculations have found that it failed for strongly correlated systems. 

mm 

The 

breakdown of the RSRG method in the strongly correlated systems is ascribed to its lack 
of account of the interactions between the blocks, i.e. merely considering isolated blocks in 
RSRG imposes wrong fixed boundary conditions for the blocks which should actually be in 
open boundary conditions. 6] 

In order to improve RSRG's performance for strongly correlated systems, various new 
schemes have been proposed. These new schemes can be sorted into two types according to 
their theoretical philosophies. One type focused on searching for new optimal criterion of 
truncation other than the energy criterion used in the traditional RSRG methods. Among 
them, the density matrix renormalization group (DMRG) method proposed by White in 



1992 
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has been shown as an extremely accurate technique in solving one-dimensional 
(ID) strongly correlated systems. The DMRG method projects the wavefunction of a 
larger block (the superblock) onto the system block and then uses the eigenvalue of the 
systems block's reduced density matrix as the basis truncation criterion. 

The other type focused on taking into account the influence of discarded excited states 
in the "isolated" block as an account of the interactions between the blocks, while the 



energy criterion of truncation is retained. 
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191 ] In order 



to achieve this goal, Lepetit and Manousakis suggested the use of a second-order quasi- 

g a new transformed 



Hamiltonian from the calculations on dimers of blocks 



ll|. To formalize the latter idea in 



a more general way, Morningstar and Weinstein introduced the contractor renormalization 
group (CORE) method in 1994. Q, Q, Q, H, Q CORE introduces effective interactions 
between the blocks from the exact spectrum of dimers or trimers of blocks by virtue of 
Bloch's effective Hamiltonian theory 20|]. In CORE, the entire system is divided into blocks 
with even number of sites and several eigenstates are kept in each block. In 2001, Malrieu 
and Guihery proposed another new improved version of RSRG - real space renormalization 
group with effective interactions (RSRG-EI) 17, Q, 19|. CORE and RSRG-EI originated 
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from similar basic ideas, but RSRG-EI uses larger blocks with odd number of sites and keeps 
only one eigenstate in each block and RSRG-EI can be applied to infinite systems due to 
the fact that it is iterative. Because CORE and RSRG-EI divide the whole system into a 
few blocks and only calculate blocks and fragments with a certain size exactly instead of 
the whole system, the increase of computational costs is only linearly scaling with the size 
of the whole system. 

Generally, despite of DMRG's great success in ID systems, the degeneracy in the re- 
duced density matrix spectrum will increase when DMRG is applied to the quasi-lD large 
systems with complicated connectivity or two-dimensional (2D) systems, which implies the 
impractical requirement of much more states to be retained to sustain the numerical pre- 
cision. In the meantime, although CORE and RSRG-EI are computationally efficient for 
large systems, the accuracy is a bottleneck due to the limitation in the fragment size that 
can be solved exactly. Meanwhile, exact or numerically high-precision calculation of larger 
blocks and fragments is actually very essential to achieve reliable account of the intra-block 
and inter-block correlation effects, and consequently a reliable account of the total system. 
In this paper, we propose a new scheme - block density matrix renormalization group with 
effective interactions (BDMRG-EI), which is a generalization of CORE with the incorpo- 
ration of DMRG. It combines the advantage of DMRG in accuracy and that of CORE in 
computational costs. 

The outline of the paper is as follows: in Section 2, the details of BDMRG-EI methodology 
are introduced; in Section 3, demonstrative computations on a ID chain and a double-layer 
polyacene oligomers are presented; finally, we conclude and summarize our results in Section 
4. 

II. METHODOLOGY 

A. Bloch's effective Hamiltonian theory 

The one-to-one correspondence between two iso dimensional subspaces can be well de- 
scribed by the theory of effective Hamiltonians established by Bloch. Firstly, one may 
consider a m-dimensional model space, onto which one would like to build an effective 
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Hamiltonian iJ'^^-^. Let us call the model space as and its projector as Pq 

Po=J2\I){I\ / = !,.. .m (1) 

leSo 

The main task of H'^-l'-f is to fulfill the requirement that its m eigenvalues are exact eigenvalues 
and its eigenvectors are the projections of exact eigenvectors in 5*0. This means that if 

H\iljk) = Ekl^pk) k = l,...n{n>m) (2) 

then H'^-ff satisfies the condition that 

H'-f^Pom =£kPom k = l,...m (3) 

The m eigenstates \ipk) which are targeted by the effective Hamiltonian span the so-called 
target space S, with the projector P, iso dimensional to 5*0. If the transform operator from 
So to S is {P = ^Pq), then the effective Hamiltonian satisfies Eq. (jHj) is 

H^f^ = PqHP = PohnhPo (4) 

On principle, one may choose S arbitrarily provide that f2 exists, but an obvious rational 
choice is to select m eigenstates with largest projections in 5*0 to span S, i.e., one should 
maximize Ylkes II l-^oV^/c) ||- After 5" is constructed, H'^-l'-l' can be obtained from the following 
equation. 

H^ff = J2\Po^Pk)£k{Poi^k\ (5) 

k 

It should be mentioned that, H^-^'-f may be non-Hermitian if |PoV'fc) loses the orthogonality 
after projection, i.e., (PoV^fcl-PoV^i) 7^ ^ki- Therefore, orthogonalization treatment of \Poipk) is 
necessary before the construction of H^^^ . 



B. DMRG 
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Over the last decade, the DMRG method [7|, l8[ has emerged as the most powerful method 
for the simulation of strongly correlated ID quantum systems. Here, we would like to give 



a brief introduction to the main physical princi_p 
be found in some recent reviews about DMRG 19 
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DMRG. More technique details can 
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1. Decimation of state spaces 



Due to the exponential growth of the number of degrees of freedom in quantum many- 
body systems, the exact simulation in the complete Hilbert space is apparently not feasible 
beyond very small sizes. Among a great number of various approximate simulation methods 
for quantum many-body systems, one important class of approximate simulation methods 
attempts a systematic choice of a subspace of the complete Hilbert space which is anticipated 
to contain the physically most relevant states. All variational and renormalization group 
techniques are within this group, and the essential question is of course to identify the best 
decimation strategy which will depend on both the system and the physical question. 

Here, let us particularly focus on the decimation strategy in strongly correlated ID quan- 
tum systems. Imagine we grow the system successively, adding site by site. The original 
system, which we refer to as an old block, is assumed to be effectively described within 
a state space la) of dimension M, the new site within a state space \a) of dimension N . 
Obviously, the state space |/?) of the new block composed of the old block and the newly 
added site will have the dimension MN , and for the prevention of exponential growth it will 
be decimated down to the dimension M. Whatever the physical decimation strategy one 
uses, the states of the new block will be a linear combination of the old states, 

\^) = = EE^"/3M|a)|a) (6) 

where matrices A of dimension M x M have been introduced, one for each |o"), such that 
the matrix elements encode the expansion coefficients: v4Q,^[cr] = {aa\f3). The introduction 
of the y4-matrices allows to encode the iterative growth of states of larger and larger blocks 
by matrix multiplications. 

Pushing the above idea further it turns out that, we can generally describe a quantum 
state of a L-site lattice emerging from a decimation procedure as 

|^)= ^M^i]^'N...A^K]ki...aL) (7) 

0"l---0-L 



Such states are referred to be matrix product states (MPS) [23|, [24, l25|, l26|, \2j\ (also known 
as finitely correlated states). Now, the main problem for simulating the state of a quantum 
system effectively becomes how we can find suitable A-matrices with limited dimensions such 
that approximates the state well. If we want to target the ground state of a Hamiltonian 
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FIG. 1: Block definition in DMRG. 



H, apparently the optimal choice for the A-matrices is to find the prescription yielding 
those A that minimize {iplHlip) with the constraint of {ip\ip) = 1. However, directly working 
out this expression leads to a highly non-linear expression for the energy in A, which is 
numerically very difficult to be solved. One way to turn the problem low-scaling consists 
in providing a starting set of A-matrices in a warm-up procedure, preferably close to the 
true solution, and then to repeat the following process iteratively: keeping all A-matrices 
in 1-0) fixed, with only one or two exception(s), and then minimizing {iplHlip) with respect 
to these one or two fiexible A-matrix (matrices) and updating these one or two A-matrix 
(matrices) according to the newly determined with minimized {iIj\H\iIj). After that one 
shifts the position(s) of the fiexible A-matrix (matrices) forth and back through the entire 
chain successively. An optimal approximation of the ground state, which is very close to 
the real global minimum of {t/j\H\ilj), can be gradually reached. Historically, this process 
corresponds to the so-called finite-system algorithm in standard DMRG language. 

2. DMRG procedures 

However, there are still some questions not answered: how can we derive A-matrices from 
the newly determined and can we retain all the important information provided by 
Now we will introduce an efficient basis truncation scheme which retains most physically 
important information and can easily derive new A-matrices. 

Assume we join one block of length n — 1 and another block with length L — n — 1 with 
two active sites explicitly inserted between these two blocks, as illustrated in Fig. [H In this 
case, the left block of length n — 1 would be described by a M-dimensional Hilbert space 
with states which can be described by A matrices as. 




(8) 
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Similarly the right block would be described by states |mf_^_^). The left block of length 
n — 1 combined with n-th site is now called new block S, and the right block of length 
L — n — 1 combined with [n + l)-th site is now called new block E. So the quantum state of 
the entire system can be described as, 

M N N M 

i^) = E E E E ^ms^4^'av^')\^^^^) 

mS = l 0-5 = 1 a^ = l m^ = l , s 

= EE^^^-i^)i^') 

i=i j=i 

where \m^(j^ = \m^a^ = \j) and = MN, = MN. During the finite-system 
DMRG sweeps, we shift the positions of the two active sites successively to update all A 
matrices for improving our wavefunction. In order to prevent the exponential growth of 
matrix dimension, one has to find a suitable truncation of the basis |z) to M states, whose 
expansion in \m^) and \a^) will just define the desired A'^[(Tn]. Similarly, we may also have 
to find a suitable truncation of the basis \j) to M states, whose expansion in \m^) and \a^) 
will just define the desired /l""'"^[cr„+i]. 

In order to find the optimal decimation strategy, one can perform Schmidt decomposition 
for li/j) = After the Schmidt decomposition, one can get the wavefunction in 

the following new form j^: 

Nsch midt 

\^)= E v^\o\o (10) 

where the scalars Wa are non-negative, and and {w^) are newly linear recombined 
vectors with Wa for S and E parts respectively. Upon tracing out the E or S part of the 
state the reduced density matrix for the S or E part is easily found to be 

^Schmidt 

ps= (11) 

a=l 

Ns chmidt 

PE= (12) 

a=l 

Apparently, Wa are the eigenvalues of the reduced density matrix for the S or E part. 

The approximate wavefunction where the space for S or E part has been truncated to be 
spanned by only M orthonormal states minimizes the distance to if one retains the M 
eigenstates of 'ps or 'pe with the largest eigenvalues ly^.j^ This is just the key truncation 
criterion in DMRG. 
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One could also look for a decimation criterion from the view of maximizing the retained 
biparticle entanglement between S and E parts under truncation. As biparticle entanglement 
is defined as = — Yla '^al-og2Wa and normally one has a large number of relatively small 
eigenvalues, this again leads to the same truncation prescription: one must retain the M 
eigenstates of 'ps or pE with the largest eigenvalues Wa- (see e.g. j^) 

For local quantities, such as energy, magnetization or density, it was also found that the 
errors are of the order of the truncation weight Cp = 1 — J2a=i '^oi^ which emerges as the key 
error estimate. 

One can therefore derive A-matrices from the effective truncation which retains only M 
eigenstates of with largest eigenvalues Wa of the reduced density matrix built by the newly 
determined \ in DMRG, the most important information for the purpose of minimizing the 
errors of the approximate wavefunction and maximizing the retained biparticle entanglement 
between the blocks is kept if one retains the eigenstates of with largest eigenvalues Wa of 
the reduced density matrix. 

Obviously, how efficiently DMRG can work depends on how quickly the ordered eigen- 
value spectrum Wa of the reduced density matrix p will decay. Empirically, in ID systems, 
density matrix spectra of gapped quantum systems exhibit roughly system-size indepen- 
dent exponential decay of Wa- So, DMRG calculations can efficiently yield reliable accuracy 
that can be comparable to exact calculations for ID quantum systems if one controls the 
truncation weight ep very small. 

C. Block density matrix renormalization group with effective interactions 

Similar to that in CORE and RSRG-EI, in BDMRG-EI, the whole system is also divided 
into a certain number {Nuock) of blocks with sites and only several low-energy states are 
retained in each block. But in BDMRG-EI, we use DMRG instead of exact diagonization 
to solve the m lowest eigenstates for each block. For example, m eigenstates with 
lowest eigen-energies e{A)i are obtained for block A with A^^ sites. 

At a later stage, we will derive the effective interactions between the neighbor blocks 
from the nearly exact spectrum of dimers or trimers of blocks according to Bloch's effective 
Hamiltonian theory. 

Let me take fragment A — B m the whole system as an example. Fragment A — B 
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contains two blocks: block A and block B, and we build small Hilbert spaces spanned by 
m lowest eigenstates for each block. Then, the model space 5*0 {dim{So) = irP) is obtained 
from the direct product of these two small spaces. The projector of .So is expressed as 
^0 = E£ ma - By^mA - B)il where - B),) = |<I>(A),) ® Then, we 

perform a standard DMRG calculation for the fragment A — B in oder to get the nearly 
exact low energy eigenstates \ip{A — B)^) with the energy e{A — B)/; for the block dimer. 
Before using Eq. to calculate the the effective Hamiltonian of fragment A — B, we should 
select largest projections on 5*0 to build the target space S and perform orthogonalization 
treatment of PolV^fc)- iWk) = J2^ M){^i\^Pk)) 

According to the above introduction to DMRG, |$(y4)j) calculated by DMRG can be 
described as 

=J2A'''{i)A'''{i) ...A^^A^i) 

(13) 

\cri) ® IcTs) (g) . . . laN^) 

To enforce the boundary conditions, one may require the leftmost matrix A^^ to be 1 x M- 
dimensional, and the rightmost matrix A^^a to be M x 1. We can write and 
liplA — B)i) similarly. 

=^S'"^A + l(z)S''^-4+2(^) . . . A'"^A + ^S (z) 

(14) 

Wna+i) ® Wna+2) ® • • • Wna+Nb) 

\ij{A - B),) = J2 ■ . . C'^A+iVfl (^) 

(15) 

(g) \a2) ® . . . laNA+Ng) 

|$(A — B)i) can be expressed as the direct product of |$(y4)j) and [^(i?)^) as the following 
equation: 

|$(^ - B),) = A"^ (j) . . . A^^A (j) ® B^^A+^ (k) ... B^^A+r^B (k) 

(16) 

Therefore, the overlap between DMRG states |$(^ — B)) and 1-0(74 — B)) can be calculated 
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easily according to the following equation: 




(C^^ . . . C^a+Nb ) 



(17) 



Ena+Nb 



where En can be calculated successively as the following formula: 



E, 



for n = 1; 



for 1 < n < Na] 



(18) 



'n 



= < 



5-"* ® En-iC' 

V 



hi NA + l<n<NA + Nb. 



for n = Na + 1; 



When the the overlap between DMRG states can be easily calculated with the above 
mentioned MPS representations, one can calculate the the effective Hamiltonian of fragment 
A — B according to Eq. @. 

Similarly, the effective Hamiltonians of other block dimers, trimers or other fragments 
can also be derived. 

After the determination of the the effective Hamiltonian of different fragments, effective 
interactions between the blocks can be determined by subtracting the contributions of all 
connected sub-clusters as the following equation. 



Finally, the effective Hamiltonian for the whole system acting on a truncated Hilbert 
space is obtained through the summation of a series of the local effective Hamiltonians and 
then directly diagonized with Lanczos or Davidson algorithm. 





(19) 




(20) 
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FIG. 2: Ground state energy of a 240 site ID chain. 



III. TEST APPLICATIONS 

As a test case of our approach, we have implemented BDMRG-EI calculations to a 240 
site ID chain and a double- layer polyacene (Pac(2)) oligomer containing 48 hexagons within 
the spin-1/2 Heisenberg model {H = JJ2{ij) '^iSj), where J is a positive exchange constant, 
and Si represents the spin operator of zth site with (ij) denoting summation restricted to 
nearest neighbors. In our test BDMRG-EI calculations, only two lowest eigenstates were 
kept for each block and only two-body inter-block actions were considered, i.e. m = 2 and 
the summation for the whole system's effective Hamiltonian in Eq. (l20l) is restricted to the 
first two terms. 



A. ID chain 

For the 240 site ID chain, we have four block definition patterns as follow: I. composing 
of 10 blocks (24 sites in each block); II. composing of 6 blocks (40 sites in each block); III. 
composing of 4 blocks (60 sites in each block); VI. composing of 3 blocks (80 sites in each 
block). 

In Fig. [21 we illustrate the ground state energies calculated by BDMRG-EI and DMRG 
with different numbers of kept block states (M). Obviously, direct diagonalization of the 
Hamiltonian matrix is impractical at this time for the 240 site ID Heisenberg chain. Mean- 
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FIG. 3: First excited state energy of a 240 site ID chain. 

while, as can be seen in Fig. [21 the DMRG calculated energy has converged with the increase 
of M. Accordingly, this converged value could serve as a reference value for the comparisons 
of BDMRG-EI results. 

As a whole, calculated energies by BDMRG-EI slow down and converge gradually with 
the increase of M, although there are some oscillations, such as the result of BDMRG(I) with 
M=400. This is due to non-variational nature of the BDMRG-EI method. But obviously, 
the range of the oscillation decreases remarkably with the increase of M and the enlargement 
of block size. 

With large M values, all BDMRG-EI calculations with different block definitions give sat- 
isfactory results in good agreement with the DMRG value, with largest error within 0.006 J. 
This implies that, for ID systems with simple connectivity, the block size in BDMRG-EI 
calculations will not affect the ground state results very much. 

In Fig. ini the calculated energies of first excited state are presented. It can be seen 
that the errors between the BDMRG-EI results with DMRG standard values for the excited 
states are much larger than those for the ground state. The values of these errors range from 
0.005 J to around 0.02 J. At the same time, the influence of block size on the computational 
accuracy are much more significant comparing with that in the case of the ground state. 
When the block size is too small, BDMRG-EI calculated energies may incorrectly fall below 
the DMRG result, and smaller block size will lead to larger deviation from the DMRG result. 
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n=8 



FIG. 4: Definitions of blocks for the BDMRG-EI calculation of a Pac(2) oligomer containing 48 
hexagons. 

This shows BDMRG-EI's non-variational nature again. Only when the block size is as large 
as composed of 60 or 80 sites, BDMRG-EI calculations can give reliable results. These facts 
elucidated that electrons are much more delocalized in the excited states, which makes it 
more difficult to evaluate the electron correlation effect in such cases. This implies that 
incorporation of DMRG into the CORE method is very necessary in calculating the excited 
states of conjugated systems. 

B. Double-layer polyacene (Pac(2)) 

The above ID chain example is a system with simplest connectivity. When the system's 
connectivity is more complicated, inter-block correlations will be more significant. As a 
consequence, the infiuence of block size on the calculation accuracy will be more remarkable. 
Here, we take a Pac(2) oligomer containing 48 hexagons as an example. In our BDMRG-EI 
calculations, three different block definitions are adopted, as shown in Fig. |H 

Calculated energies of the ground state and the first excited state of this Pac(2) oligomer 
are shown in Fig. [5] and Fig. [6l Similar to the situation in ID chain, DMRG calculated 
values have gradually converged with the increase of M for this Pac(2) oligomer, and the 
DMRG result with M = 500 serves as the reference value for comparisons. 

In BDMRG-EI calculations with small block sizes, the errors are around 0.2J for both 
the ground state and the first excited state, which are definitely intolerable. However, we 
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FIG. 5: Ground state energy of a Pac(2) oligomer containing 48 hexagons. 




FIG. 6: First excited state energy of a Pac(2) oligomer containing 48 hexagons. 

can still find from the figures that the results will improve remarkably with the increase of 
the block size. As an example, the errors in BDMRG-El(Vll) for the ground state and the 
first excited state fall within 0.05J. These values are acceptable for such complex system. 

It should be mentioned that, even the smallest block size (6-8 hexagons in BDMRG- 
EI(V)) in BDMRG-EI calculations are far beyond the applicabilities of the normal CORE 
or RSRG-EI methods(at most 2-3 hexagons). 



14 



IV. CONCLUSIONS 



This work proposes BDMRG-EI method as an improvement of the traditional CORE. 
Within the basic framework of CORE, BDMRG-EI incorporates DMRG calculations into 
the solution of the blocks and block dimers or trimers. 

Our demonstrative calculations showed that the block size effect in systems with compli- 
cated connectivity is much more significant than that in systems with simple connectivity, 
and that in the excited states is also much more significant than that in the ground state. 
Since the DMRG calculations can achieve high accuracy for systems much larger than those 
could be exactly treated, BDMRG-EI could take much more intra- and inter-block corre- 
lations into account comparing to the original CORE and RSRG-EI methods. Hence, the 
results of CORE are greatly improved by BDMRG-EI. Demonstrative calculations of a 240 
site ID chain and a Pac(2) oligomer containing 48 hexagons convinced that high accuracy 
can be achieved in BDMRG-EI calculations when enough large block size is chosen. 

Another advantage of BDMRG-EI is its efficiency with low computational costs. Because 
it solves only the blocks and fragments with a certain size by DMRG instead of the whole 
system, BDMRG-EI provides a new low-scaling electron structure method which can be 
applicable to very large systems with simple or complicated connectivity. 
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